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Abstract 

Recently an analytical model was presented that treats the nonlinear signal distortion from the Kerr nonlinearity in optical 
transmission systems as additive white Gaussian noise. This important model predicts the impact of the Kerr nonlinearity in 
systems operating at a high symbol rate and where the accumulated dispersion at the receiver is large. Starting from the suggested 
model for the propagating signal, we here give an independent and different calculation of the main result. The analysis is based 
on the Manakov equation with attenuation included and a complete and detailed derivation is given using a perturbation analysis. 
As in the case with the published model, in addition to assuming that the input signal can be written on a specific form, two 
further assumptions are necessary; the nonlinearity is weak and the signal-noise interaction is neglected. The result is then found 
without any further approximations. 

Index Terms 

Optical fiber communication, communication system nonlinearities, nonlinear optics, wavelength division multiplexing. 
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I. Introduction 
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Q ■ TN AN OPTICAL transmission system operating at a high symbol rate, the propagating signal is rapidly evolving due to 
^ | _1_ chromatic dispersion (CD). It has been found numerically that after a relatively short distance of propagation, the probability 
O ■ density function of all four quadratures of a polarization- multiplexed signal become Gaussian with zero mean and a variance 
5a ' related to the signal power JT]. This can also be shown analytically by using the central limit theorem and is intuitively 
£>~v understandable since the dispersed signal at every point in time can be viewed as a coherent superposition of many signal 
pulses. During the propagation, a nonlinear phase shift is induced in proportion to the local power by the Kerr nonlinearity. 
As the CD is compensated in the receiver, this phase shift will give rise to a residual signal distortion. Numerically, it has 
been observed that this distortion is very similar to additive white Gaussian noise (AWGN) [1|. This observation is of great 
importance as it suggests that if no attempt is made to compensate for the nonlinear effects, then modeling the nonlinear signal 
distortion as AWGN is possible. 
. In 2011, an analytical model was presented that calculates the noise-like nonlinear distortion for a quite general wavelength 
y—t ' division multiplexing (WDM) system Q, J5). This signal distortion was named nonlinear interference (NLI) and recently an 

■ extensive paper on the same topic was published [4|. The work reported here is closely related to these publications but it 
' should be pointed out that there are many results in the literature that address the same or similar questions. To name one 
. example, the work for OFDM by Chen and Shieh [5 1 has many similarities in both the approach and the results. However, the 
' previous work in the area seems to be well described in the introduction of J4] and we will not further elaborate this topic 

r j~^ here. 

J> ' In the publications |2]]-|[4]], a suggestion is given for how to model the signal, the four-wave mixing (FWM) of the different 
K^l ' signal spectral components is calculated, and the corresponding power spectral density (PSD) is found. However, the derivation 
' of the main result is quite short. Partly this is because known FWM results are used flSJ. In this paper, an independent detailed 

■ derivation of the resulting PSD is carried out. The calculation starts from the Manakov equation with power gain and attenuation 
included and we use the signal model suggested in [2|. The calculation is based on a perturbation approach previously used, 
e.g., to investigate intrachannel cross-phase modulation (XPM) and intrachannel FWM, which gives rise to "ghost pulses" in 
systems using on-off keying |7|-[9|. This approach has later been used both to analytically study systems, see for example iflOl . 
Ifm . and to compensate for the NLI, see for example Ifl2l . Using the perturbation approach, we here present a self-contained 
calculation of the PSD of the NLI directly from the model equation. This work also goes beyond [4| since, as described in 
Section III-B1 we perform the calculation for a more general system. 

The organization of this paper is as follows: In Section [II] the perturbation analysis is introduced and a formal solution that 
is valid for all input signals is given. In Section [III] we find the solution corresponding to the specific input signal suggested 
in (2, which is then used in Section [TV] to calculate the NLI, i.e., the PSD of the perturbation. We then calculate the PSD 
corresponding to the result in [0] Eq. (18)] by studying a specific system choice in Section [V] Finally, we conclude. 



Pontus Johannisson is with the Photonics Laboratory, Department of Microtechnology and Nanoscience, Chalmers University of Technology, SE-412 96 
Gbteborg, Sweden, e-mail: pontus.johannisson@chalmers.se. 
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II. Perturbation analysis 

The calculation of the NLI in assumes that the FWM is a weak effect. This assumption is stated explicitly as "the 
pump is undepleted". We proceed in a similar way by introducing the complex envelope of the electric field in the x and y 
polarizations according to A = (A x ,A y ) T and writing this as A = A; + A p , where A; is the linearly propagating signal, 
i.e., the signal in the absence of any Kerr nonlinearity, and A p is a small perturbation. This implies that the nonlinear effects 
may not become significant and this constitutes the first assumption of the model. The second assumption is that the input 
signal can be written on the form suggested in J3]. The range of validity for this signal model is a separate question and we 
do not discuss this question here. A third assumption is that the signal-noise interaction is neglected as the calculation of the 
NLI does not involve amplifier noise in any way. This is also an assumption of [2| and in order to account for amplifier noise 
under this assumption, noise corresponding to the total AWGN from the amplification is added just before the receiver. These 
three assumptions are sufficient and no further approximations are necessary to carry out the analysis. 



A. The perturbation equation from the Manakov equation 

In order to describe transmission using polarization multiplexing, we start from the Manakov equation |[T3l and include 
power gain and attenuation. It should be noted that the Manakov equation is obtained by averaging over the polarization 
rotations which are assumed to be fast and this equation does not take polarization mode dispersion into account. Denoting 
the group-velocity dispersion by fa(z), the power gain by g(z), and the power attenuation by a(z) we have 

*f - - 7W(A H A)A + ^A, (1) 

where A H A = \A X \ 2 + \A y \ 2 is the sum of the power in the x and y polarizations and the nonlinear parameter 7(2) = 
(8/9)(fco^2/^-eff)- In the expression for 7, fco is the wavenumber corresponding to the center frequency, 112 is the Kerr 
coefficient, and yl e ff is the effective area of the optical fiber. The power gain g(z), which is set up using erbium-doped fiber 
amplifiers (EDFAs) and/or Raman amplification, is assumed to have no frequency dependence, i.e., the gain is flat over the 
bandwidth of the signal. The perturbation is introduced according to 

A x ,i \ / A x ,p 

where Ai(z,t) solves the linear equation obtained by setting 7 = in (fl}. We find that 



A <+ A » = (7 + IT'). w 



(A H A)A = (1/1,., + A T ,„f + \A yl + A,X) ( A :-\ + A /-' I . (3) 



Ay^l ~t~ Ay,p 

Our intention is to study A x and we will consider the two cases that either (i) \A x j | and \A y j | are of the same order of magnitude 
(transmission using polarization multiplexing) or (ii) A y j = (single-polarization transmission). The first assumption above 
can then be strictly formulated as l^^pl <C \A x j\ and \A y ^ v \ <C |A<;,z|- This allows us to approximate the nonlinear term to 
leading order according to 

(A H A)A«(|^| 2 + |^, ; | 2 )f^V (4) 



Inserting this into the Manakov equation, we obtain 
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( A x,l " 


\~ A X p\ 


_ fa 


d 2 


I A x j + A^.pN 


dz 


\ A v,i ~ 


\~ A y,pJ 


2 


dt 2 


\Ay,l + A y,pJ 



y,l 



, 1 1 1 1 \A X {\ ,g — a (A x , 1 + A xp > 
l(\A x ,i\ + \A V: i\ )[*•)+ [ a , a P ), (5) 



y AyJ J 2 \Ay.l ~\~ Ayp 

where the z-dependence is now implicit for compactness. The fact that A; solves ([TJ when 7 = 0, implies that 

■ d (A x .p\ fad 2 (A XtP \ |2 . |2 . (A x A . g-a ( A x . p \ 

[ A ) =^^72 A -7(kWI +\Ay,i\ ) , ' A • 



'dz \Ay.pJ 2 dt 2 \A y ,pJ ,VI x ' 11 1 v ' n J \AyjJ 2 \A y , p 
We will study the x-polarized perturbation, which is described by 

8A X r, fad 2 A xv ., . . , .,. . q — a , 

~ihT = ^t~d^ + + 1^ 1 )A *> 1 + -^r A ^v- (7) 

For compactness, we temporarily introduce S(z, t) — A x j \ 2 +\ A y j | 2 ) A Xi i to denote the source term in the partial differential 
equation for the perturbation, which is then written 

dA XtP .fa d 2 A XtP g-a 
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B. Formal solution 

We now derive a general formal solution without any assumptions about the input signal or the system parameters. To 
describe the power evolution, we introduce P(z) as a function that satisfies the equation 

dP 



dz 



[g[z) - a(z)]P. 



(9) 



In the absence of any Raman amplification this function varies as e~ az between the amplifiers. The EDFA gain can be modeled 
by a (5-function in g(z) to obtain the discontinuities in P(z) at each z corresponding to the location of an amplifier. Introducing 

A x , P {z,t) = \J P{z) ip(z, t), we obtain the equation 



—VPip + Vp^ + i y Vp^t - — -Sp* = s, 



or 



+ 1 



d 2 i> s 



dz 2 dt 2 yfp 
We proceed by describing the accumulated dispersion, B(z), by a function that satisfies 

dB 



dz 



= P 2 (z). 



(10) 



(11) 



(12) 



If there is lumped dispersion compensation (such as a chirped fiber Bragg grating), then B will have discontinuities. This can 
be modeled by a <5-function in /^(z) to obtain 

(13) 



(14) 



B{z)= / p 2 (QdC 
Jo 

Fourier transforming^ the equation, and denoting this operation by tilde, we obtain 

... n2 fo r S 



dz 



VP' 



which can be written 



Integrating f Q ■ d(, we use 



to obtain 



dz 



(^-i{2nf) 2 B/2^ _ e -i(2TTf) 2 B/2 



S 



P 



-i{2-nf) 2 B/2 



-1 <= z 



1> 



C=o 



= {V»(o,/) = o} = 



_ n\ - „-t(27vf) 2 B( Z )/2 



$( Z , /) = e ^ffB(z)/2 f Z e ^(2,f) 2 B(0/2S{^fl 

Jo VW) 



(15) 



(16) 



(17) 



The final perturbation at the receiver is described by ip(L, f), where L is the total system length (possibly containing many 
spans with fibers and amplifiers). Assuming that perfect EDC is performed, the exponential before the integration is canceled 
and we obtain 



c -i(2irf) 2 B(z)/2 S ( Z J) j y 



(18) 



Assuming that the power at z = L is equal to the power at the transmitter, Pq, we have A x , p (L,t) = y/P(L) ip(L, t) — 
\fPo~il){L,t) and we can write the general expression for the perturbation as 

'^-^BM^SM^^ (19) 



where p(z) — P{z)/Pq has been introduced. This function describes the normalized power evolution through the system. 
Using the definition of S, we find 

A x ,p{L>, f) - 



L p -i(2nf) 2 B(z)/2 







-T[i^(\A X:l \ 2 + \A Vtl \ 2 )A x , l } dz, 



(20) 



where T denotes Fourier transformation. 

'We use the same definition for the Fourier transform as Proakis [14j, i.e., u(f) = f^L, u(t)e~ t2w ^ t dt and u(t) = u(f)e l27T f t df . 
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III. Signal model 

Provided that the first assumption is fulfilled, (l20t is the solution to (HJ without further approximations. However, we need 
to choose the system parameters, select the boundary condition at z = 0, i.e., the input signal, and perform the integration. 
When trying to do this with an accurate modeling of the signal pulses, the calculations become cumbersome and the final 
expression needs to be averaged over the data. We here instead use the second assumption, i.e., we write the initial field in 
the way suggested in |2|. This model is 

oo 

A x ,i(0J) = Vh E WG*(fc/oW-fc/o), (21) 

k= — oo 



A/,i(0,/) = V^ E CkJG v (kf )S(f-kfo), (22) 



k= — oo 



where and (k are complex independent Gaussian random variables of unit variance and the input signal PSD of the x and y 
polarizations are denoted by G x (f) and G y (f), respectively. In the following, we will suppress the infinite summation limits 
for notational convenience. We account for dispersion and power variations during the propagation by modifying this to 

A x ,i(z,f) = vME^VW)^/ - kf )e^ k M 2B ^'\ (23) 



A y ,i(zJ) = \ffoW) E C^G y (kf ) S(J - fcio)e l(2 ^ /o)2s(z)/2 . (24) 
fe 

Using this linear solution, we can now calculate the perturbation using (l20l >. In order to do this we need the expression 

K,lAy 



FM\A Xil \ 2 + \A ytl \ 2 )A x ^ = i 1 T[Al l A* l ] + i^A^AyM. (25) 



The solution ( f2Qb can be written 

A xp (L,f) = i [% e -i^I)'B/2 ^M dz + l f\ e -^ ffB/ 2 ^ A ^ A U d (26) 

Jo VP Jo s/P 

where we again omit the z dependence for notational compactness. We have 



A x ,i(z,t) = E & \^ G M) e i <- 2 " cfo W" k M Bl \ (27) 



Ay,i{z,t) = y/faJ2W G v( k M e *™ t0)teii2 * k,0) > ™ 



and get 



Al,A* x , = ( y/f^^Hy/G^e^f^e^^j & VC^Io) e^K^f^j 

x fv^ECv/G^^e-^^^e-^ 2 ^) 25 /^ , (29) 



A x ,iA y ,tA* yil = v^E^^^) el(2 ^ /o)tel(2,rfc/o) 72 ^pY.^ G y^)^ lSo)t ^ lh 



,fB/2 



x ^Ecv^M) e " t2 ™ /o)fe " 1(2 ™ /o) " B/2 j- ( 3 °) 

These expressions can be simplified to 

AhA*, = (/op) 3 / 2 E &6CVG,(fc/o)G,(;/o)G,Wo)e !2 ^ +i - m)/ ° i e l(2 ^ (,:2+,2 - m2)/ » B/2 , (31) 

k,l,m 



A x ,iAy t iA* tl = (/ p) 3/2 E a(iCv / G#/o)G J (i/o)G J (™/o)e a(H ^ )/ ' ! e ,(2 *) 2(t2+|2 - m2)/ « B/2 (32) 

k.l.m 



The Fourier transforms are 

FlAl^] = (/op) 3 / 2 tei&y/G x (kf )G x (lf )G x (mf )6(f -(k + l- m )h)e^^^- m ^ B l\ (33) 

k,l,m 

F[A Xtl A y>l A; tl ] = (/op) 3 / 2 E (kOC^G x (kf )G y (lf )G v (mfo)S(f -(k + l- m )f Q )e^ 2(k2+l2 - m2 ^ B / 2 . (34) 

k,l,m 

The integrands of ( |26i > can then be written 

ie -i^ f fBi ^\ A lxK^ = ie ^ f?B /2 f m p £ a6CVG,(fc/o)G x (// )G x (m/o) 

x <5(/ - (k + I - m)f a y^ 2 ( k2+l2 - m2 VZ B l\ (35) 
ie ^ rfB/2 HA x ,A = 7e - j(2OT/) ^ /2/ 3/ 2p ^ aoC^( fc /o)^(^o)^(m/ ) 

x *(/ - (fc + / - m)/ ) e l(27r)2(fe2+/2 - m2)/ ° B/2 . (36) 
We can move the exponential functions containing B after the summation signs to obtain 

f\A2 A* ] 

7e -W)-B/ a L ^ = /o 3/2 7p 5- W C>/G,(*/D)G.(!/o)G,(m/o) 

x 5(f — (k + l — m )f Q y(^f(k 2 +i 2 -^)fZB/2 e ^(2^(k + i- m ffZB/2^ (37) 



V p fc,/,m 



We use that 



to obtain 



x 5(i -(k + i- m )f ) e i(^r(k^- m ^B/2 e -i^r(k+i- m r^B/2_ (38) 

(k 2 + l 2 - m 2 ) -(k + l-mf = -2(k - m)(l - m) (39) 

7e -i W fB/2 ^5j = /o 3/2 7p £ e^CVGx(fc/o)G a a/ )G s (m/o) 
v p fc,i >m 

X(J(/-(fc + I" m )/ ) e -^ 2 (fc-m)((-m)/ 2 S (4Q) 

7e-^^^/ 2 ^^^ = / 3/2 7P E MC^x(fc/o)G,(!/o)G ( ( m / ) 

xS(f-(k + l- m)f a )e- lA " 2{k - m ^ l - m ^'o B . (41) 



The solution can therefore be written 



P (L, f) = i [ fo /2 JP E ZkZi&y/G x (kf )G x (lf )G x (mf )6(f -(k + l- m )/ ) e -^ 2 Cc-™)(*-™)/ 2 i3 dz 

+ * / L /o /2 7P E \/ G - ( fe /o)G, (Z/ )G V (m/ ) <S(/ - (* + Z - m )/ ) e -^ 2 ( fc -™)(<- m )/o B dz. (42) 



Inverse Fourier transformation gives 

.3/2 



J ° k,2,m 

+ */ 3/2 / L ip E aocv /G ^ fc /°) G ^^ o ) G ^ m / o ) e42 " (fc+ ^ m)/ote_ ^ 2( ^ m)(i " m)/o2s dz - < 43 > 
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We can also rearrange integration and summation to obtain 

A XtP (L,t) = if* /2 ]T 66C\/Gx(fc/o)G,(i/o)Gx(m/o)e i2l ( H, - m » Ai f 7pe -^ 2 ^-m)(i-m)fSB dz 



k.L: 



*/o /2 E ^OCV / G :c (fc/o)G y (Z/o)G y (m/o)e^( fc+i -™)^ t / 7pe -^ a (fc-m) (I -m)/ 2 B ^ (44) 

fej.m 



Considering j(z), p(z), B(z), and L to be given system parameters, we introduce 

C klm = C(kf , lf , m/ ) = / lpe -^ 2 (^)(i-m)f B dz 

Jo 

to write this as 

A x , p (L,t) = ifj 2 J2 ^ m 6aCv / Gx(yo)G I (!/o)G I (m/o)e & ( t + ! - m »/» f 



(45) 



k.L 

,3/2 



z/ 3/2 E &™aOCV^( fc /o) G 9 ('/o)G 9 Wo)e* 



(k+l-m)f t 



k.L: 



k,l,m 

For later use we notice that 



= z/ 3/2 £ C fcime t2 " (fc+ '-" l)/o *6 VG,(fc/ ) y/G x {lf Q )G x {mf ) + OC^(I/o)G„(m/ ) ) . (46) 



7pe 



-i4Tr 2 (k-m)(l-m)fgB 



dz = jpdz^'ypL, 



(47) 



where p and 7 are the maximum values of p(z) and 7(2:) in the interval z £ [0, L}. This means that |Cfcz m | is upper bounded 
by a constant as k, I, m, and /o are changed. 

IV. Power spectral density 

When the input signal is modeled by (l2Tl i and d22l >, then the perturbation is given by d46i >. The next step is to calculate the 
corresponding PSD. Using the Wiener-Khinchin theorem lfT4l p. 67], we have that the PSD is the Fourier transform of the 
autocorrelation according to 



-OO 



where 



(48) 



(49) 



R = E{A x , p (L,t 1 )A* XiP (L,t 2 )} 

and E{-} denotes the expectation operator. As will be shown, the latter expression can be written as a function of the time 
difference r = t\ — £2- We will suppress the infinite limits for notational convenience. We have 



A XtP {L,t x )A*{L,t 2 ) 



(50) 



*/o /2 E Cki m e t2 < k+l - m ^t 1 ^ Gx{kfo) U&y/GMWzimfo) + 0Cn^G y (lh)G y (mf ) 

k,l,m * 

-z/o /2 E Cw m e- a< * +l '- mVo Hl-y/G x {k%) (cU'VG x (l%)G x (m%) + CC m ^G y (l'f )G y (m'f ] 



/o 3 E C fclm C l V m ,e ,2 ^ +1 -")^V- 2 "( t ' +l '- m ^» i2 aSVG I (yo)v/G I (^/o)x 



k A ,rn 



ti&y/G x (lfo)G x (mf ) + QC m ^G y (lf )G y (mf )) UUm'V G x (l%)G x (m%) + Q( m ^G y (l'f )G y (m%) 



We see that the complete autocorrelation function will consist of four terms 

R = Ri + i?2 + R3 + R4 



(51) 
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where 



\ k,l,m 



(fc+!-m)/ ti e -»27r(fc'+J'-m')/o*a 



fe'J',) 



x ^iC m C k ^U^G x {kf )G x (lf )G x {mf )G x (k'f )G x (l'f )G x (m%) , (52) 



i? 2 = E \f* C klm Cl, Vm ,e^ k+l - m ^e 



{k+l — m) fati —i2ir(k'+l , — m r ) f Q t 2 



k,l,m 



x ^ie m e k <i< m ^G x (kf )G x (lf )G x (mf )G x (k%)G y (l'f )G y (m'/o) j , (53) 

D _ TP J f3 p n* A2TT{k+l-rn)f Q t 1 -i27v{k'+l'-rn')f Q t2 

^3 WO / ^ ^klm^k'l'm'^ c 

\ k,l,m 
k',l',m' 

x efeOCS^emV G -( fc /o)^(«/o)G y (m/o)G K (fc7o)G4«yo)Gx(m'/o) 1 , (54) 



z2jr(fe+(—ro)/o*l -i2ir(k'+l'-m r )fot 2 



x ^CiCMCm^^^^ (55) 

These four terms in the autocorrelation will give rise to four different terms in the total PSD, which we denote by G XiPj i, 
G XlPt 2, G x>Pt 3, and G XiPt 4, respectively. We need to find these four terms individually. 

A. The random variables 

It is seen that in (T52b— (T55b. all terms except the £ and £ are deterministic. In order to carry out the expectation operation, 
in this section we therefore need to investigate these random variables. 

The £ and £ are complex independent Gaussian random variables of unit variance, i.e., E{£fc} = 0, E{£ 2 } = 0, E{|£fc| 2 } = 1, 
Vfc. Analogous expressions hold for Using this, we can simplify the six-dimensional sums (|52T>— (T55T> in the following way. 
First assume that one of the summation variables, say k, has a value different from all other summation variables, i.e., k is 
unique. Due to the independence, we then have 

Ei^CmCk^U'} = e{&> m&m>u>} = o. (56) 

An identical argument holds if some of the £ are replaced by £. This implies that the expected value is always zero when 
one of the summation variables is unique. Thus, we have the condition that no summation variable can be unique. Second we 
assume that no summation variable is unique, but there are three pairwise equal values, where each pair has a unique value. 
Assume for example that k = I, k' = V, m = m'. Then 

Eu fc 6c^*U'} = mkuuuum} = e{&&} mut} nc m u} = ml} nm 2 } n\u\ 2 } = o. (5?) 

We conclude that each random variable must be paired up with the complex conjugated version of the same random variable. 
This conclusion reduces the dimensionality of the summation to three or less. However, the dimensionality cannot be less than 
three, because then the expression G x p — > as f — > 0. To see this we first notice that 

\G x ,p\ < |G XjPi i| + |G X)Pl 2| + |G XiPi 3| + |Ga;,p,4|, (58) 

and since all terms behave similarly in this respect, we can study, say, G x , P ,i- Let us select the one-dimensional case k = I = 
m = k' = I' = m'. We then have 

i? a =Elf^\Ckkk\ 2 e^ khT \^k\ 6 GUkfo)\ =/ 3 ^|C fcfcfc |V 2 ^ T E{|^| 6 }G3(fc/ ). (59) 

I k ) k 
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We get 

G x , pA = f^J2\ C kkk\ 2 n\Ck\ 6 }GUkf )6(f - kf ). (60) 

k 

We see that G XiP; i — > as / — > by identifying the Riemann sum and writing the expression as 

G x , P ,i = fl J \C(h,h,h)\ 2 mk\ 6 }Gl(h)S(f h) dh = .fo 2 |C(/, /, f)\ 2 n\tk\ 6 }Gl(f). (61) 

Remembering that the value of |C| 2 is upper bounded by a constant, the fact that G XiPi i — > is now obvious. An analogous 
argument can be made for the case of a two-dimensional sum. We conclude that the summation must have dimension exactly 
three, i.e., each random variable must be paired up with the complex conjugated version of the same random variable and all 
three pairs must have unique values. 



B. The first term in the PSD 

We now study the first expression, i.e., G XiPi \(f) = F[R\{t)\. Following the above rules for how the indices can be chosen, 
we have six possible combinations: 



(62) 



k 


= k' 


1 = 1' 


m = to' 


k 


= k' 


I = to 


m'= I', 


k 


= V 


I = k' 


TO = TO' 


k 


= V 


I = m 


m'= k', 


k 


= m 


I = k' 


m'=V, 


k 


= m 


1 = 1' 


m'= k'. 



We need to study these different possibilities individually. 

1) The case k = k', I = V, m = to': We then have 



{ k,l,m 

= Jo 3 E \Cki m \ 2 e l2 < k+l -^ foT G x (kf )G x (lfo)G x (mfo), (63) 



G x , P ,i(f) - fo E \Cki m \ 2 G x (kf a )G x (lfo)G x (mfo)S(f - (k + I - m)/ ). (64) 

k,l,m 

Letting fo —¥ 0, we get 

G x , pA (f) =111 \C(fi, / 2 , h)\ 2 G x {.h)G x {h)G x {hW - h - h + h) dhdhdh 

|C(/i, / 2) fi + h - f)\ 2 G x {h)G x {h)G x {h + h - /) dhdf 2 . (65) 



2) The other cases: Compared to the case above, the case k = V, I = k', m = m' is obtained by swapping k' and /'. Since 
Ckim = Cikm and the rest of the expression is clearly invariant under this change, this yields an identical result as the case 
above. The other four cases are all equivalent to the case k = k', I = m, to' = I'. We get 

I HV J 
I HV J 

= fl E C m Ct n ,e^ kf '' T G x (kf Q )G x (lf )G x (l%), (66) 



kW 



G x , P Af) = flY, C ^ C ^ G -( k h) G -( l ^)G x {l'h)5{f-kf ). 

Ml' 



(67) 
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Letting f — >• 0, we get 

G x , P ,i(f) = jjj / 2 , / 2 )C*(/i, / 3 , h)G x {h)G x {f 2 )G x {h)6{f - h) dhdhdh 

= jf C(f, f 2 , f 2 )C*(f, f 3 , h)G x {f)G x {f 2 )G x {h) df 2 dh. (68) 

However 

C(/i, / 2 , / 2 ) = C(/i, / 3 , /s) = C(/ 2 ,/i, / 2 ) = C(/ 3 , /i, / 3 ) = / 7(«)pW dz = Co, (69) 

Jo 

which is a real number. This gives 

G x , pA (f) = C%G x (f) J I G x (f 2 )G x (f 3 ) df 2 dh = ClG x {.f) J G x (f 2 ) df 2 J G x (f 3 ) df 3 = C%P*G x (f). (70) 

3) The total PSD for the first term: The six possible index selection cases split into two groups of degeneracy two and 
four, respectively. We get the total expression 

G x , P ,i(f) = 2 J J |C(/i, f 2 , h + h - f)\ 2 G x {h)G x {h)G x {h + h - f) dfidf 2 + AClP 2 x G x {f). (71) 

C. The second term in the PSD 

We now study the second expression, i.e., G Xi p^(f) = F[R 2 {t)}. Following the above rules for how the indices can be 
chosen, we have two combinations 

k = k' I = m ml =V 
k = m I = k' m'= I' 



(72) 



The reason that we have fewer possibilities is that the expression contains both £ and (, which are independent. 

1) The case k = k', I = m, m'= I': We then have 

[ k,l,V 

= / 3 J2 CmCl w e^ kiaT G x {kf )G x {lf )G y {l'h), (73) 

k,l,V 

G x , p , 2 (f) = fo E CkiiC* m ,G x (kf Q )G x (lfo)G v (l%)8(f - kf ). (74) 
k,i,v 

Letting / — > 0, we get 

G x , p , 2 (f) = HI C(h,f 2 ,f 2 )C*(f u fs, h)G x {h)G x {f 2 )G y {h)5{f - h) dhdf 2 df 3 
= C 2 jj G x (f)G x (f 2 )G v (f 3 ) df 2 dh 

= C 2 P x P y G x (f). (75) 

2) The case k = m, I = k', ml '= I': We notice that this case is obtained from the above case by swapping k and I and will 
therefore give the same result. 

3) The total PSD for the second term: We get the final expression 

G XtPt2 (f)=2CiP x P y G x (f). (76) 

D. The third term in the PSD 

We now study the third expression, i.e., G x , p ^(f) = ^[^(t)]. Following the above rules for how the indices can be chosen, 
we now have two combinations 

k = k' I = m m! =V 

k = V l = m m'=k' { ' 

The first case is identical to the first case for the second term in the PSD. The second case is obtained from the first case by 
swapping k! and V and will therefore give the same result. We get 

G,,,,,:/: -1C:: I >,!),( (78) 
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E. The fourth term in the PSD 

We now study the fourth expression, i.e., G XiP! a(J) = J-[Ra{t)]. Following the above rules for how the indices can be 



chosen, we have two combinations 



k = k' 1 = 1' m = ml 
k = k! I = m ml ' = I' 



k.Ll' 



k.LV 



Letting f — > 0, we get 

G x , pA (f) = 1 1 1 C{h,hJ 2 )C*{h,h,h)G x (h)G y {h)G y {h)5{f - h)dhdf 2 dh 



(79) 



1) The case k = k', 1 = 1', m = m': We then have 

R±=e\ /o 3 ]T l^™| 2 e^ (fc+i - m)/o Ha| 2 |0| 2 |C m | 2 G,(fc/o)G,(?/ )G,(m/o) 

i k,l,m 

= / 3 E \Ckim\ 2 e l2 * {k+l - m)foT G x (kf )G y (lfo)G y (mf ), (80) 

k,l,m 

G x , pA (f) =foJ2 \Cki m \ 2 G x (kf )G v (lfo)G v (mfo)S(f -(k + l- m)/ ). (81) 

k,l,m 

Letting f — > 0, we get 

G x , pA (f) = /// \C(fi, / a , h)\ 2 G x {h)G y {h)G y {h)5{f - h - fr + f 3 ) dhdf 2 df 3 

\C(fi, /a, fi + /a - f)\ 2 G x (h)G v (f 2 )G v (h + / 2 - /) dhdf 2 . (82) 
2j 77ie case k = k', I = m, m' = I': We get 



[ k,i,v 

= / 3 ^ C m C* klll ,e^ k ^G x (kf )G y (lf a )G y (l'f ), (83) 



G XlJ) , 4 (/) = / 3 E c mC* m G x {kh)G y {lf )G y {l'h)5U - kf ). (84) 



C(/, /a, /a)C*(/, / 3j h)G x {f)G y {f 2 )G y {h) df 2 df 3 
= ClP*G x {f). (85) 
5) 77ie total PSD for the fourth term: We get the final expression 



G X , P M) = JJ \C(fi, /a, A + /a - fWG x {h)G y {h)G y {h + / 2 - /) rf/id/a + C l ^G x {f). (86) 
77ie foffl/ PSD 

We can now find the total PSD by summing up all the terms according to 

G XtP = G XtPi i + G x , p , 2 + G XiPt 3 + G XjPt 4 

= 2 ff \C(f l7 f 2 , h + h- })\ 2 G x (h)G x {f 2 )G x {h + /a - /) dhdf 2 



+ JJ \C(A,h, h + h~ f)\ 2 G x (h)G y (h)G y (h + h-f) dfidh 

+ Cl{^P 2 x + 4P x P y + P*)G x (f). (87) 

We notice that the x polarization acting on itself is twice as effective as the y polarization acting on the x polarization. As is 
clear from the derivation, the reason for this is the level of degeneracy. This is similar to the case of the double influence of 
XPM as compared to SPM, when the nonlinear Schrodinger equation is approximated by a coupled system with one equation 
for each WDM channel [15 p. 264]. The final term of (l87l i is of no consequence for the transmission. The origin of this term 
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is the phase modulation from the entire propagating field acting on itself. In practice, we expect this to just cause a rotation 
of the received signal constellations. However, in the perturbation analysis, this type of phase modulation gives rise to these 
extra terms in the PSD. As (|87| i can be considered the main result of this calculation, we repeat the expression for C for ease 
of reference 

/a, h + h -!)= t i(z)p{z)e'^ 2 ^^-^ dz. (88) 
Jo 

It should be noticed that C is determined when the physical parameters of the channel are selected and \C\ 2 is a measure of 
the FWM efficiency. Then, by choosing the input signal PSD we obtain the PSD of the NLI. 



V. Transmission system example 

We now calculate the resulting expression for a system of the type considered in [2 |. The system consists of A>smf spans, 
each containing a standard single-mode fiber (SMF) followed by an EDFA. There are no dispersion-compensating fibers (DCF) 
or any other optical dispersion compensation. Instead, the CD will be compensated for by using DSP in the receiver. The first 
span starts at z = zq = and the last span ends at z = zn smf . The fiber parameters a, /3 2 , and 7 are assumed to have no 
z-dependence. Furthermore, we assume that identical signals are launched in the two polarizations, i.e., we set G y = G x and 
we remove the terms that are not due to FWM to rewrite (1871 as 



Gx,p 



|C(/i,/ 2 ,A + /2 - f)\ 2 G x (h)G x {f 2 )G x {h + h - f)dhdf 2 . 



(89) 



We have B(z) — 02 z and temporarily introducing k — 47r 2 (/i — f)(f2 — /) we find 

/a, fx + /a - /) = / iv^ B[z) dz 



JO 

i^i SMF 



-a(z-z„_i) — iKfiiz 



dz 



E 

71=1 
-/VsMF 



-a(z n 



-') 



Lsmf „ — iK/3 2 "isMF 



7 



-(1 



1 



-ZK/3 2 (n-l)L S MF 



-iK ( 32-LsMF-'VsMF 



where we also used the assumption that all SMFs have the same length, denoted by Lsmf = z n — z n -\. We find 



(90) 



|C(/i,/a,/i + /a-/)| 2 =7 2 



a + inf3 2 



1 



3 — JK^2-£'SMF-'VsHF 



g — aLsMFg — »K&isMF 



a + 1KP2 



sin 2 (k/3 2 ^smf^V S mf/2) 
sin 2 (K/? 2 isMF/2) 



(91) 



We then find 



G x ,p = 37 



-aL SM F-i477 2 (/ 1 -/)(/ 2 -/)/3 2 L s 



a + z47r 2 (A-/)(/ 2 -/)/3 2 
x G x {h)G x (f 2 )G x {h + /a - /) rf/irf/a- 



sin 2 [2 7 r 2 (/ 1 - /)(/ 2 - /)/3 2 L S MF^SMFj 



sin 2 ^ 2 ^-/)^-/)^^] 



(92) 



In order to compare this expression with the coherent expression in [4 Eq. (18)], we must account for the fact that a different 
convention for the Manakov equation is used in [4]. Thus, we need to replace 7 with 87/9. Furthermore, the NLI PSD in [4. 
Eq. (18)] is the summation over both polarizations, i.e., we here get Gnli = G XtP + G VtP = 2G XtP . Similarly, we here have 
Gtx = G x + G y = 2G X . (From [4, Eq. (13)], it seems like Gj x is equal to G x . However, this is not the intended meaning of 
Gtx- We thank the authors of |4) for explaining this.) Finally, there is a difference in how the attenuation is defined. For the 
definition used here, we refer to ((T) and ((9). When these differences in notation are taken into account, the results become 
identical. 
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VI. Conclusion 

We have presented a derivation of the power spectral density of the nonlinear interference under the assumptions that 
(i) the nonlinear effects are weak, (ii) the signal can be written as suggested in J2), and (iii) the signal-noise interaction can 
be neglected. Using these three assumptions, we obtain the PSD expression for a general system. Applying the result to a 
system consisting entirely of SMF and accounting for the differences in notation, we obtain a result identical to that presented 
in @-B|. 
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